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Abstract. We analyse a period spanning 35 years of activity in the Sao Paulo 
Stock Exchange Index (IBOVESPA) and show that the Heston model with stochastic 
volatility is capable of explaining price fluctuations for time scales ranging from 5 
minutes to 100 days with a single set of parameters. We also show that the Heston 
model is inconsistent with the observed behavior of the volatility autocorrelation 
function. We deal with the latter inconsistency by introducing a slow time scale to the 
model. The fact that the price dynamics in a period of 35 years of macroeconomical 
unrest may be modeled by the same stochastic process is evidence for a general 
underlying microscopic market dynamics. 
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1. Introduction 

In the last decades much attention has been devoted to the modeling of asset returns 
by the Finance community, the major drive behind this effort being the improvement 
of pricing techniques for derivative contracts Ij. The pricing problem is amenable 
to analytical solution for some stochastic volatility models |2] such as Hull- White [Hj, 
Stein-Stein and Heston |5j. Despite differences in methods and emphasis, the cross 
fecundation between Economics and Physics, which dates back to the early nineteenth 
century (see [H] and [7]), has intensified recently [Sj. Following the tradition of statistical 
physics, substantial effort has been made to find models capable of elucidating the 
basic mechanisms behind recurrent features of financial time series such as: returns 
aggregation (probability distributions at any time scale) CHI, volatility clustering 
[TT] . leverage effect (correlation between returns and volatilities) [121 CHI, conditional 
correlations [TH HHj . 

Recently, a semi-analytical solution for the Fokker-Planck equation describing the 
distribution of log-returns in the Heston model has been proposed jU]. The authors 
were able to show a good agreement between the return distributions of a number of 
developed market stock indices and the model, for time scales spanning a wide interval 
ranging from 1 to 100 days. 

In this paper we show strong evidence that the Heston model is capable of explaining 
the return distribution of the Sao Paulo Stock Exchange Index (IBOVESPA) in a period 
that span 35 years of political and economical unrest, with hyperinflation periods, 
currency crises and major regulatory changes. We also show that the Heston model 
can explain the diffusive process of IBOVESPA from minutes to months with a single 
set o parameters. However, we observe that the Heston model is inconsistent with the 
measured volatility autocorrelation function and we introduce an extension to the model 
along the lines discussed in which exhibits the correct behavior. 

2. The Heston Model 

The Heston Model describes the dynamics of stock prices St as a geometric Brownian 
motion with volatility given by a Cox-IngersoU-Ross (or Feller) mean-reverting 
dynamics. In the Ito differential form the model reads: 

dSt = Stfitdt + St^tdWoit) (1) 
dvt = - 7 [vt - 9]dt + K^/n[ dWi{t), 

where Vt is the volatility and dWj are Wiener processes with 

{dW.it)) = 0, {dWjit) dWkii)) = [5jk + (1 - 6,k)p] 5{t - t) dt. (2) 

The term ^Jvl avoids unphysical negative volatilities, 6 is the macroeconomic long term 
volatility and /it represents a drift also at macroeconomic scales. 
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Interestingly, v{t) = Yl'^=i -^ji^) ^ Feller process (P) if Xj are Ornstein-Uhlenbeck 
processes (OU) defined as [T7]: 

dX,it) = -^X,it)dt + ^dW,it), (3) 

where dWj describe d independent Wiener processes. This observation restricts the 
possible parameter values to ci > 2 as the probability density of the volatility must 
vanish at zero to be consistent with the empirical data. To see how the volatility 
process in (P) emerges from OU processes we apply Ito's Lemma to get: 

d d 2 

dvt = -bdtJ2^j + Yl (4) 

j j j 

Using the definition of v and the properties of the Wiener processes it follows that: 
'd 



dvt 



— bvt 



dt + a^tdW. (5) 



The volatility process in (0) can be recovered by a few variable choices: a = n, b = 'j 
and e = f — . Note that, given the dimension d, there are only two free parameters in 
(P. The OU processes in Q may be regarded as the primary microscopic sources of 
volatility and the condition of non- vanishing volatility implies that a = ^ > 1. 



3. The Fokker-PIanck Equation Solution 

We now outline the solution of the Fokker-PIanck equation (FP) describing the 
distribution of log-returns proposed by Dragulescu and Yakovenko P . 

As we are mainly concerned with price fiuctuations, we simplify equation ^ by 
introducing log-returns in a window t as r{t) = ln(S'(t)) — ln(S'(0)). Using Ito's lemma 
and changing variables by making x{t) = r{t) — dr^ir we obtain a detrended version 
of the return dynamics that reads: 

dx = -^dt + ^tdW^. (6) 

The FP equation is, therefore, given by: 

dP d Id d"^ 1 d'^ 

that has to be solved for the following boundary condition P{x, v, \ Vi) = 6{x)6{v — Vi). 

A Fourier transform in x followed by a Laplace transform in v leads to a 
partial differential equation of the first degree that can be solved by the method of 
characteristics [T^ : 

d 



d 



Q = ~iPvQ, (8) 



dpv 

where F = ipi^Px + 7 and the boundary condition is Q{px,Pv, \ Vi) = e"^"^* 
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The unconditional distribution of log-returns can be obtained by inverting the 
Laplace and Fourier transforms and integrating first over the volatility v and then over 
the initial volatility Vi. 



Pfix) 



dvi P*{vi) 



dv 



27r 



(9) 



'0 JO 

where P* is the stationary solution for the volatility distribution (see |;9 for details). 
Integrating © we finally get: 



Pix) 



e exp 



a ( y - In cosh - ^t{j>x, 



where 



O = a/F^ + k2 {pI - ip^ 
^t{Px) = In 



fi2_r2 + 27F , /Qt 

1 + — tanh — 

27^ V 2 



(10) 

(11) 
(12) 



4. Volatility Autocorrelation and Leverage Functions 

The formal integral of (0) is given by: 

vt = {vq - 6) e-^i* + 6 + ^1 dWi{u) e-^(*-") 

Jo 

A simple calculation gives the stationary autocorrelation function: 



C(r I 7, 6, k) = lim 



{VtVt+r) - {Vt){Vt+r) 



a 



The leverage function is [T!^ : 



Tf I n N _ r {dxt{dxt+rf) 
L{t 7, e, K, p) = hm 



PkH{t)G{t) e-^\ 



where dxt is given by is the Heaviside step function and: 



G'(- 



Vt exp 



(13) 



(14) 



(15) 



(16) 



To simplify the numerical calculations we employ throughout this paper the zeroth order 
appoximation G{t) ^ G{0) = 9^^. The approximation error increases with the time lag 
T but is not critical to our conclusions. 

In the next section we simultaneously fit (UH), (fT^ and the model ffTUI) to 
fluctuations of the BOVESPA index at a wide range of time scales. 



5. Fitting IBOVESPA data 
5.1. The data 

Three data sets were used. IBl consists of daily data from IBOVESPA inception on 
January, 1968 to December, 2002. IBl was adjusted to take into account eleven divisions 
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Figure 1. Left: Large deviations in daily log-returns with identifiable exogenous 
causes: New York Stock Market Crash (1); Economic Plans (2-4); Presidential crises 
(5), Recovery from Mexican, Asian and Russian crises (6); Asian and Russian crises 
(7) and Currency crises (8). Right: (Top) Macroeconomic trend in the log-prices. 
(Bottom) Trend and fluctuations in the last ten years. 



by 10 introduced in the period for disclosure purposes pUj. IB2 consists of daily data 
from January, 2001 to August, 2003. IBS consists of high-frequency data from March 
3, 2001 to February 14, 2003 and from June 6, 2003 to August 26, 2003. 

5.2. Trimming 

Large deviations are explained both by endogenous dynamics and by exogenous 
shocks. An adequate microscopic model for the price formation process has to explain 
spontaneous large movements and the system response to external shocks. However, here 
we only intend to model the typical behavior at mesoscopic scales (price dynamics). In 
emerging markets, large movements with exogenous identifiable causes are frequent. In 
order to obtain a reliable fitting for the typical fluctuations we expunged from the data 
set large deviations connected to major structural changes, in Figure 1 we identify those 
externalities. For a brief historical account of recent Brazilian economy see |19j . 

5.3. Filtering the Macroeconomic Dynamics 

In Equation ini we extract the macroeconomic drift to focus on price fluctuations. This 
drift at long time scales reflects the effect of inflation, economic growth, business cycles 
and the riskless interest rate. To extract we have employed a low-pass Savitzky-Golay 
smoothing filter of degree two |21| with averaging over a four years moving window {i.e. 
one Brazilian presidential term). In Figure 1 we show dr fir for the whole data set 
(top) and the fluctuations around the trend in the last ten years (bottom). This choice 
of smoothing technique is heuristic and may introduce artifacts at time scales of lengths 
comparable to the moving windows (4 years in our case). We, therefore, restrict our 
analysis to a maximum period of 100 days and leave a thorough analysis of this aspect 
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Figure 2. Top: Volatility autocorrelation: IBl (circles), Heston model (dashed) and 
two time scales (solid line). Bottom: Leverage function: IBl (circles) and approximate 
Heston model (see text, solid line). 



to another occasion. 



5.4- Fitting Parameters 

After filtering the trend we have to fit four parameters: the long term mean volatility 
6', the relaxation time for mean reversion I/7, the volatility fluctuation scale /t and 
the correlation between price and volatility p. It became apparent in ^U] that a simple 
least squares fit that takes into account only (|T0|) yields parameters that are not uniquely 
defined. The task of finding parameters fitting all the main stilyzed facts simultaneously 
is not trivial. In this paper we do not intend to focus on parameter estimation, hence 
we adopt a heuristic procedure and do not provide rigorous error bars. These questions 
will be addressed somewhere else. We first estimate the long term mean volatility 9 
directly from the daily log-returns as: 

TV 



e 



N 



(17) 



where x*-^-* stands for daily detrended log-returns. 

A proper fit must be also consistent with the non- vanishing volatilities constraint, 
f > ^ > 1. We, therefore, minimize a constraint free cost function and adjust k to 
attain consistency with the constraint. 

The following cost function is employed: 

2 



1 



L{Ti I 7,61, K,p) 



+ 



1 ^ 

1=1 



(18) 
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IBOVESPA Daily 68-02 (IB1) 




X (log-return) 

Figure 3. Probability densities of IBOVESPA returns from 1 to 100 days for the 
dataset IBl. The vertical scale is logarithmic and is multiplied by a constant for 
better visualization. 



The empirical leverage function is [13): 



where x 



(1) 



In 



-Y 



M 



t=l -^t 



(1) /^(l) 



X 



t+T 



02 

fit and the empirical autocorrelation function is: 



(19) 



M ^ 



t=l 



X 



(1) 



X 



(1) 

t+T 



(20) 



We then employ the parameters obtained above to fit the probability density of returns 
and adjust the parameter p to the empirical data minimizing the mean squared error. 
We will discuss a more systematic approach to parameter estimation elsewhere as our 
main aim here is restricted to showing the description capability of the Heston model. 
The parameters found are in the following table: 



parameter 


IBl 


IB2-IB3 


e 


7.8 X 10"^ days'^ 


5.2 X 10"^ days-^ 


1/7 


9.0 days 


5.8 days 


P 


-0.20 


-0.15 


K 


1.3 X 10^2 days-^ 


1.1 X 10-2 days^^ 


d 


2.03 


2.98 



Figure 121 shows IBl data and curve fits. It is clear that, under the vanishing volatility 
constraint, a single relaxation time is not simultaneously consistent with leverage and 
autocorrelation functions. The best fit for the autocorrelation is shown as a dashed 
line. We deal with this inconsistency in the final sections of this paper by extending the 
Heston model to include two relaxation times. 
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Figure 4. Probability densities of returns at 5, 10 and 15 minutes, overnight and lunch 
time jumps (dataset IBS). 



5.5. Intraday Fluctuations 

At first glance, it is not clear whether intraday and daily returns can be described by the 
same stochastic dynamics. Even less clear is whether aggregation from intraday to daily 
returns can be described by the same parameters. To verify this latter possibility we 
have to transform units by determining the effective duration in minutes of a business 
day Teff. The dimensional parameters are, therefore, 6'*^^^^ = O/Tefj, 7*^^^^ = ^jT^ff 
and K^^^^ = n/Teff. In order to avoid the effects of non-stationarity we also use data 
from the same period (IB2 and IBS). From IB2 and IBS we found the effective duration 
of a day to be about Tg// = 540 minutes. The Sao Paulo Stock Market opens daily at 11 
a.m. and closes at 7:S0 p.m. local time, resulting in exactly 510 minutes. The effective 
duration infered from data exceeds, therefore, in SO minutes the real duration of a normal 
business day, an amount that is explained by the effect of overnight and lunchtime jumps. 
In Figure 4 we show a histogram of jumps compared to intraday returns at 5, 10 and 15 
minutes. Considering both overnight and lunchtime jumps, between 10 to SO minutes 
must be added to a normal business day, just the right amount to make T^ff compatible 
to the hypotheses of a single set of parameters describing time scales from minutes to 
months. In Figure El we show observations and data for IB2 and IBS with parameters 
given by the table above. These findings are suported by similar results presented in 
|22] for high-frequency returns of individual stocks traded in NASDAQ and NYSE. 

The datasets IBl and IBS yield T^ff = 600 minutes, the 90 minutes excess is 
explained by the well known non-stationarity of the parameters revealed in the table 
above, that shows that the mean variance in the S5 years period (IBl, 6 = 7.8 x 10~^ 
days~^) is higher than the mean variance observed in the last two years (IB2, 9 = 
5.2 X 10~^ days~^). It is interesting to stress that in an emerging market like the 
Brazilian structural changes related to political events are frequent. In Figure IHl we try 
to put the the macroeconomic non-stationarity of the mean variance ^ in a historical 
perspective, identifying some important political events. The robust statistical behavior 
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IBOVESPA High-Frequency (IB3) IBOVESPA Daily 01-03 (IB2) 




X (log-return) x (log-return) 



Figure 5. Left: Probability densities for intraday returns from 1 minute to 4 hours. 
Circles are empirical data from IBS while lines correspond to the theory. Right: 
Returns from 1 to 100 days. Circles are empirical data from IB2. A single set of 
parameters (see Table) is used through all time scales considered. The vertical scale is 
logarithmic and is multiplied by a constant for better visualization. 



here reported are particularly surprising when one considers the extremely abrupt 
structural changes that take place in an emerging market. 

6. Multiple Time Scales 

It is clear in Figure El that the Heston model is not capable of describing the 
autocorrelation function. To rectify this inconsistency we propose the inclusion of a 
second, slower, time scale into the dynamical model along the lines of . We assume a 
stochastic dynamical model where the stochastic volatility reverts to a second stochastic 
volatility with much longer relaxation time 71 ^72- The new model reads: 

dS = S{t)i2dt + S{t)y^ dWo{t) (21) 
dv = --fi[v{t) - e]dt + Ki^/vi{} dWi{t), 
dd = -j2[0{t)-9o]dt + K2^/9it) dW2{t), 

where 72/71 << 1 e dWj are Wiener processes defined as: 

{dWM = 0, (22) 

{dWj{t) dWk{i)) = CjkS{t - i) dt, 

where C is a correlation matrix with Cjj = 1, C12 = p and C13 = C23 = 0. 
In this model the autocorrelation function acquires the following form: 

C{r) = + , (23) 

«! a2 

with a2 = where 6 stands for the average of 6 given 60 . In Figure El we show a 
fit of this autocorrelation function to the data, the new relaxation time is 7^^ = 144.9 
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Figure 6. Historical perspective for the non-stationarity of the mean volatihty 
calculated using the dataset IBl and a two years wide moving window. The numbers 
identify the main political events as follows: military dictatorship ends after two 
decades (1), first presidential election (2), presidential impeachment (3), second 
election (4), reelection (5), third election (6). 

days and K2 = 1.0 x 10~^ days~^. It is possible to solve the Fokker-Planck equation for 
this extended model in the limit where 02 ^ ai- We will discuss this model in detail 
elsewhere |22|- 

7. Conclusions 

We showed that the Heston model can reproduce the diffusion process of a 35 years long 
time series of IBOVESPA returns at a wide range of time scales from days to months. We 
also show that the Heston model can explain the aggregation of returns from minutes 
to months with a single set o parameters that change in the macroeconomic scales. 
However, the Heston model, with a single relaxation time, is not capable of explaining 
the behavior of the autocorrelation function. We, therefore, proposed an extended 
version of the model with the addition of a slow stochastic dynamics that yields an 
autocorrelation function consistent with the observations. 

It is surprising that a single non-trivial stochastic model may be capable of 
explaining the long-term statistical behavior of both developed and emerging markets, 
despite the known instability and high susceptibility to externalities of the latter. 
We believe that this robust non-trivial statistical behavior is evidence for more basic 
mechanisms acting in the market microstructure. Perhaps the search for underlying 
common mechanisms (or laws) that can explain empirical data is the main contribution 
of Physics to Economics. This contribution might be particularly useful to the field 
of Econometrics in which a common view is that a theory built from data 'should be 
evaluated in terms of the quality of the decisions that are made based on the theory' 
Clearly, these two approaches should not be considered as mutually exclusive. 
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